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DESCRIPTION OF RESEARCH 

Our study of the time variability of the galactic X-ray sources 
which exhibit nonperiodic intensity fluctuations has concentrated on the 
source Cygnus X— 1. The data which we have examined were obtained with the 
wide f ield-of-view collimator aboard the Uhuru satellite. Our research pro- 
gram, however, has had three principal goals: first of all, to develop the 

computer software and expertise to perform analyses of the time variability 
of the galactic X-ray sources; second, the specific application of these 
programs to the study of the source Cygnus X-l in order to obtain a further 
insight to the emission mechanism of this source; and finally, to develop a 
better understanding of what may be learned from further studies along the 
same lines. We feel that we have succeeded in all three endeavors. 

I. Software Development 

Computer software has been written to accomplish the following: 

(a) To detrend the Uhuru data; i.e., to remove the angular response 
of the collimator from the raw data. These programs are, in fact, applica- 
ble to data received by any scanning receiver . 

(b) To calculate autocorrelation functions. Autocorrelation func- 
tions are calculated, for the detrended data for each independent observa- 
tion of an X-ray source. The energy bandwidth can be selected. 

(c) To calculate a cross-correlation function. The data in two 
operator-selected energy bandwidths are cross correlated for each observa- 


tion of an X-ray source. 
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(d) To perform Fourier analysis. The autocorrelation functions can 
be multiplied by a variety of standard spectral windows (Bartlett, Parzen, 
and Tukey) , and a power density spectrum is computed. 

(e) To provide averages. Software to determine the average overall 
independent observations of the source for any of the above functions. 

(f) To perform long-term time variability analyses. This software 
is used to calculate the unsmoothed, sample power spectrum of the long— term 
(days) variations of any chosen parameter such as the mean intensity per ob- 
servation, the cross-correlation at zero time delay, etc. 

(g) To simulate with Monte Carlo techniques sources exhibiting differ- 
ent categories of time variability and to produce probability distributions 

of parameters derived from the raw data. One of the most important aspects 
of data analysis lies in the understanding of the statistical uncertainties 
produced by chance fluctuations in the experimental data. This is particularly 
complicated in the analysis of the nonperiodic fluctuating X— ray sources where 
one is faced with two categories of signal fluctuations, viz., that due to the 
source itself and that due to photon statistics. Clearly errors in derived 
quantities are not simply found by the propagation of the estimated uncertainty 
of the latter. Furthermore, we are dealing with a situation where the prob- 
ability density function governing the source variability is not known <x pv'L— 

It has been extremely useful for us to simulate the results of experi— 
tients assuming a particular source variability such as white— noise and shot — 
noise sources. .In addition to allowing us to derive probability distribu- 
tions, this software was used to check the analysis routines for those cases 
where the results are well understood, e.g., the autocorrelation function for 
a DC sourde with counting statistics should be (on the average) zero for all 
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time delays other than zero and should exhibit fluctuations with variance 
l/N where N is the total number of time delays. Finally, this software 
allowed us to study the impact on our results of possible systematic effects 
such as, for example, the assumption that the angular response of the colli- 
mator was nontriangular , etc. 

JJ. Cygnus X-l Results 
A. Short-Term Time Variability 

Seventy-one independent observations of the X-ray source Cygnus X-l, 
obtained with the wide field-of-view proportional counter aboard the Uhuru 
satellite were examined. A detailed discussion of our analysis of the short 
term (<15 sec) variability is contained in a paper, "Short-Term Time Vari- 
ability of Cygnus X-l," which has been accepted for publication in the 
Astrophyslcal Journal (Letters) . A preprint is included as Appendix I to 
this report. 

Figure la in the Appendix shows the X-ray intensity as a function of 
time for a typical observation. All the data were obtained with time reso- 
lution of 0.192 sec and the energy bandwidth of observation selected by us 
for analysis was from 2.1 to 16.4 keV. Briefly, our results are as follows: 
(1) The time averaged autocorrelation function add its corresponding power 
density spectrum (the latter, apart from the effects of smoothing, is the 
Fourier transform of the former) are clearly not consistent with a steady 
(DC), white noise, nor periodically pulsed source. (2) Furthermore, the 
shape of the autocorrelation function, essentially a simple exponential 
with e folding time of 0.45 sec is a classic example of that which would 
be produced by a randomly pulsed or "shot-noise" source. This apparent be- 
havior would be a natural consequence of the formation of local "hot" spots 
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in an accretion disk, which most theoretical models of this source invoke. 
(3) We find no evidence for any energy dependence in the autocorrelation 
function or power— density spectrum. In fact, we have discovered a serious 
error in the interpretation of previous results by other workers in this 
field which led them to present, incorrectly, evidence for an energy depend- 
ent effect. (4) We find a 100% correlation in time on the scales of the 
order of 0.2 sec between flux in different energy bandwidths. This result 
would indicate that the hot spots or flares that give rise to low- and high- 
energy photons are identical and have the same spectrum. (5) When we com- 
pare our results with other similar analyses obtained since the spectral 
transition of 1971 March, once the error discussed in our paper is taken in- 
to account, we find that the results are in agreement and would therefore 
indicate that there is a steady feature of this source, namely, the mean 
rate of occurrence of flares and the mean cooling time. These results are 
consistent with the qualitative features of most of the standard accretion- 
disk models which invoke X-ray emission from two distinct regions: an opti- 
cally thin, inner region from which the high-energy flux emanates, and an 
optically thick, outer region from which the lower-energy flux is emitted. 
The X-ray flux in the 2— 16-keV range that we have examined is clearly from 
a single region as the 100% cross-correlation indicates. This result is 
consistent with the hypothesis that during the spectral transition the rel- 
ative intensities from the two regions shifted dramatically. The similarity 
between the posttransition spectrum and the high-energy component of the 
pretransition spectrum would further indicate that we are observing the flux 
from the inner region. The cooling time that we detect can then be used to 
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estimate the electron density of the emission volume. We wish to point out 
that if the hypothesis of two distinct emission regions is correct, then we 
would not expect to observe a 100% cross correlation between the low- and 
high-energy components of the pre transition data. We discuss this point in 
more detail in a proposal for continuation of support being submitted under 
separate cover. 

B. Long-Term Time Variability 

In addition to examining the short-term time variations of Cygnus X-l , 
we have also examined the long-term (~days) time variations of several param- 
eters over the 12 days spanned by the 71 observations. Of especial interest, 
of course, are variations occurring at the 5.6-day period detected from the 
optical counterpart. Unfortunately, the limited data span of 12 days and the 
poor quality of the Uhuru observations (in the sense that counting statistics 
dominate each individual observation) make it difficult to draw any definite 
conclusions. Nevertheless we have proceeded to examine the long-term time 
variations of the following variables. 

(1) The Background. 

The sample spectrum of the background for each observation is shown in 
Figure 1. This spectrum is perfectly consistent with white noise and indi- 
cated to us that, despite the nonuniform data sampling, we were not subject 
to gross systematic effects in the resulting power spectrum. This conclusion 
was also confirmed using the Monte Carlo simulations with white-noise sources. 

(2) The Mean Count Rate per Observation 

The mean count rate per observation for these observations, corrected 
for the spacecraft attitude, is shown in Figure 2. 
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FIG. 1. — Power spectrum of the background for 71 observations of Cygnus X-l in the time interval 

a-v 

from 9 to 22 January 1972. 
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Even in this rather limited data set we see variations in flux by 
as much as a factor of 4. The three error bars shown in the figure repre- 
sent, respectively, the typical statistical uncertainty based on the flux 
variations in each observation (cr^), \the typical uncertainty in the intensity 
resulting from the uncertainty in the aspect correction (c ? ), and finally, the 
uncertainty based on the scatter of the data about the mean (a ) . The fact 
that the latter is almost three times larger than the former two indicates 
of course that these intensity variations are not consistent with those expect- 
ed from a DC source with counting statistics. This can be seen more clearly 
from the resulting power spectrum in Figure 3 where there are three peaks, 
each of probability of chance occurrence of 5 x 10~ 3 or less. In fact, these 
results make extremely unlikely that the long-term variations follow a white- 
noise process (with variance o ) as the probability of obtaining such peaks 
by chance is inordinately small. We are then faced with the following possi- 
ble conclusions: (a) the long-term variations are consistent with those due 

to a white— noise source, and we happened to have observed the highly improb- 
able chance occurrence of three statistically significant peaks in the power 
spectrum, or (b) in general, the source is a white-noise source and two or 
even all three of the peaks in the power spectrum are indicative of a peri- 
odic process. Unfortunately, we have no a -priori reason for singling out 
any of the periods of 3.6, 2.8, and 1,6 days for any special significance. 

We do find it interesting, however, that the 2.8-day variation, shown super- 
imposed on the data in Figure 4, is in phase with the 2.8-day variations ob- 
served in the intermediate band b-magnitude of the visible companion HDE 226868 
(as reported, e.g., by Lester et al . , Nature Phys . Sci. 241 , 125 [1973]). 
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FIG. 4. — Intensity as a function of time. The 2.8-day variability is shown as a solid litie. S 
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(c) The most reasonable interpretation, however, is that the long-term 
variations are describable by some unidentified stochastic process. In 
fact, one could speculate about a "superflare 11 process as the high power 
at low frequencies could indicate. 

We should note that it is quite possible that the uncertainties in 
the aspect correction are significantly underestimated. These data were 
obtained after the Uhuru star sensor became inoperative, and spacecraft 
aspect was obtained by comparing the intensity between back— to-back obser- 
vation of the two detectors with differing fields of view. The short-term 
shot-noise variations of the source, though small compared with the errors 
produced by counting statistics, were not taken into account in the aspect 
corrections sent to us by the group at the Smithsonian Astrophysieal Obser- 
vatory. 

(3) The Shot Parameters 

Clearly a very interesting set of variables is the rate of occurrence 
of the shots (A) and their decay time T. Unfortunately, the signal-to-noise 
ratio of the Uhuru observations is quite poor so that it is effectively im- 
possible to make a meaningful measurement of these quantities for a single 
observation of the source. Thus we are forced to lump observations together 
which, of course, implies a reduced sensitivity to long-term variations. 

The resulting mean shot-model parameters as a function of phase for. 
four phases are shown in Figure 5. There is clearly a correlation between 
the value of the parameters and binary phase outside the statistical varia- 
tions, but there are not enough data to warrant an unambiguous selection of 
this period other than on a pri-ori, grounds . 
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III, Summary and Conclusions 

We have prepared and tested computer software necessary for the 
analysis of the time variability of X-ray sources. We have examined in 
detail 71 Uhuru observations of the black-hole candidate Cygnus X-l. We 
have shown that the X-ray emission can be best represented as arising from 
randomly occurring X-ray flares with a mean rate of occurrence of 18 flares 
sec ~~ and a mean cooling time of approximately 1/2 sec. Furthermore the 
data show that photons in the 2— 16-keV bandwidth arise from within the same 
flares. Finally, we have found evidence that the mean flare rate and cool- 
ing time have been constant since the source underwent a spectral transition 
in 1971 March. 

We have examined the long-term variations of several parameters asso- 
ciated with the X-ray emission from this object with the basic result that 
further detailed investigations are necessary. 
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SHORT-TERM TIME VARIABILITY OR CYGNUS X-l 
M. C. Weisskopf and S. M. Kahn 
Columbia Astrophysics Laboratory, Columbia University 

and 

P. G. Sutherland 
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ABSTRACT 

Autocorrelation functions, power density spectra, and 
cross-correlation functions were calculated for 71 observations 
of Cyg X-l obtained by the Uhuru satellite. The average auto- 
correlation function and power density spectrum were found to 
be independent of energy and consistent with a shot-noise source 
with a characteristic time constant of 0.45 s. The data in dif- 
ferent energy bands were found to be 100 percent correlated. 

Subject headings : 

I. OBSERVATIONS AND DETRENDING 

In this Letter we present the results of a study of the short-term 
time variability of the source Cyg X-l. The data were obtained with the wide 

O 

field of view (5.1) proportional counter aboard the Uhuru satellite (Giacconi 
et at. 1971). Figure la shows the X-ray intensity as a function of time for 
data obtained in the energy bandwidth of 2.1 to 16.4 keV from a single obser- 
vation of the source on 1972 January 9. 



2 


In order to remove the angular response of the detector, the data 
were examined visually to determine when the source was in or out of the 
field of view. The off-source data were then averaged to determine the back- 
ground level for the observation. After subtracting background, the on-source 
data were fit by means of a (weighted) least-squares technique to a triangular 
response function with the width, height, and center of the triangle left as 
free parameters. The weighting factors were based on the errors in each 
data point due to counting statistics, including the effects of background 
Subtraction. The detrending process was completed by subsequently dividing 
the source data by the (normalized) best-fit,, triangular response function. 

The error in each data point attributed to counting statistics was propagated 
through this piocedure. The detrended data for the 1972 January 9 observation 
are shown in figure lb. As indicated in the figure, only that portion of the 
detrended data for which the propagated counting statistics error is within 
a factor of two of the minimum value was then selected for further analysis . 

Seventy-one observations of Cyg X-l obtained by the Uhuru satellite ' 
in the time interval from 1972 January 9 27 were detrended in the manner 
described above. 


II. THE AUTOCORRELATION-FUNCTION- ANALYSIS 
In order to characterize the short-term time variability of the source, 
an autocorrelation function (ACF) p(u) was calculated for each observation where 

p(u) = r(u)/r(o) (X) 

and 


r(u) = 


N - u (S - S) (S . - S) 

Y i i + u 

L rr 

i = 1 i 


a. 

i + u 


N 

I 

i = 1 


l/o\ 


a. 

i 


( 2 ) 
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Here and are the (detrended) counting rate and corresponding error in 
til 

the i time bin, and 


N 


S = 


l s jo: 

i = 1 1 1 

N 

I l/a.‘ 

X = 1 L 


(3) 


is the mean counting rate per 0.192-s bin. 

The weighting factors CJ_^ were introduced into equations (2) and (3) 
to account for the systematic increase in the errors outwards from the center 
of the triangle. 

The average ACF for the 71 observations and its corresponding power 
* 


density spectrum are shown as the dashed lines in figure 2. It is important 
to note that the functions shown as the dashed lines in figure 2 are not 
representative of the behavior of the source but are biased due to the effects 
of counting statistics. 

More explicitly, the ACF defined in equation (1) is an estimate of the 
statistical function 


p(T) = 


<z(t) z(t + t)> _ <ZZ> T 

V "V 

z z 


(4) 


where the stochastic variable z is the detected counting rate with total 


The power density is defined as 


L - 1 

P(f) = 2A[1 + 2 l p(k) w(k) cos 2TTf(k)] 
k = 1 

where 0 < f (2A) \ A is the sampling time interval (0.192 s) , and w(k) is 
a smoothing function. We note that P(f) summed over all (N) independent, 
positive frequencies is NA. 
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variance V and the average indicated by the brackets is per'formed over an 
ensemble of observations (see, for example, Jenkins and Watts 1968). If the 
signal z(t) is stationary (in the statistical sense), then p is only a function 
of T , as indicated, and may be estimated through equation (1). Due to the dis- 
crete nature of the photon signal, we can write z = x + y with x the "true" 
source signal and y representing the effects of counting statistics. Then 
we have 


P(T) 


<zz> t <xx>^ + <xy + yx> T + <yy>^ 
__ - - 


(5) 


z z 

As the counting statistics are uncorrelated with the source, then in the 

ensemble average the cross terms vanish, and V = V + V . Since y is ef- 

z x y 

fectively a white noise signal, <yy>^ = V but <yy>^ = 0 for all x r 0, and 
thus 


p(T) - 


<xx> 


V 


V 

x r \ 

V P x (T) 
z 


with p (t) the ACF of the source. 


T - 0 
T i 0 


( 6 ) 


Equation (6) clearly shows that the measured ACF will be biased by 
the ratio of the source variance to the total variance. Only in the limit 
where the total and source variances are identical, i.e., when the fluctu- 
ations in the data produced by counting statistics are negligible, will the 
ACF defined in equation (1) represent the true time variability of the source. 
This point seems to have been overlooked in a previous analysis of six Uhuru 
observations of Cyg X-l by Brinkman et at . (1974) and led to erroneous con- 
clusions as discussed below. 
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The procedure that we adopted to determine the ACF representative of 
the behavior of the source is as follows: The total variance for each obser- 

vation is given by r(o) (see eq. [1]), and the contribution to the total 
variance due to counting statistics was estimated to be 

(7) 

The source variance is then given by - V - V . We then averaged the 71 
ACF 1 s and corrected for the effects of counting statistics by multiplying 
p(x t 0) by V (V ) where V and V are the total and source variances 

Z X Z X 

averaged over the 71 observations. We have tested this approach by means of 
Monte Carlo shot-noise simulations of the experimental data, and we find : that 
by this technique we are able: a) to estimate the source variance to better 

than 20 percent and b) to reconstruct the correct ACF with similar accuracy. 
The correction factor for the experimental data is 1.6, The corrected ACF 
and corresponding power density spectrum are shown as the solid lines in 
figure 2. The error bars are not based on counting statistics but rather 
are computed from the scatter about the mean of the appropriate variable 
for each observation. The increased errors in the corrected functions follow 
from the uncertainties in the average variances. 

A. Discussion 

The ACF shown in figure 2 falls exponentially to zero with an e-folding 
time of 0.45 s and clearly confirms the now well-established result that the 
short-term time variations of Cyg X-l are consistent with neither a steady, a 
white noise, nor a pulsed source. 

Power density spectra were determined by smoothing the ACF with a Tukey 
window w(k) = 0.5[1 + cos (7Tk/L)], where L specifies the truncation point. We 


V = N/ 


N 

I iAq' 

l 1 = i 
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have also computed power density spectra utilizing Bartlett and Parzen windows 
and find no significant differences among the results. The power density spec- 
trum shown in figure 2 was determined with L = 24, which seemed to realize a 
reasonable compromise between bandwidth and variance for these data. 

The "peak" in the power spectrum at low frequencies should be inter- . 
preted only as an indication of high power at low frequencies and not an 
indication of the presence of a pulsed component. The slight turnover near 
zero frequency is a result of the fact that P(o) - 0 and that with L = 24 
we are averaging over only a few frequencies. As the truncation point' L is 
reduced, the turnover disappears as expected. 

As was first pointed out by Terrell (1972), ACF’s and power spectra 
of the type shown in figure 2 are examples of those which would be produced 
by a randomly pulsed or "shot noise" source. This apparent behavior would 
be a natural consequence of the formation of local "hot" spots (flares) in 
an accretion disk (see for examples Thorne and Price 1975) which have a 
lifetime of the order of a second. Regardless of the physical interpretation, 
a shot process can be used to summarize the results. The data that we have 
examined are best represented by an exponential shot model with characteristic 
time of 0.45 s and a shot rate of 18 shots s ^ . 

B. Energy Dependence 

We have also examined the average ACF as a function of energy. In 
particular, we have compared the ACF for data in the 2.1 to 5.2 keV band 
with that for data in the 5,2 to 16.4 keV band. The measured ACF, equa- 
tion (1), leads to the results reported by Brinkman et ai. (1974), namely, 
that the high-energy ACF falls more rapidly to zero than that of the low- 
energy data. However, as different amounts of flux were detected in the 
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two bandwidths , one might well expect the ratio V /V to be different. In 
this case, conclusions drawn from the measured ACF’s may be invalid. For the 
data that we have analyzed, the detected high-energy intensity was 37 percent 
of the low-energy flux, and the ACF correction factor V /V was 2.8 for the 
high-energy data as compared to 2.0 for the low. Once the corrections to the 
ACF's are applied as shown in figure 3, we find no evidence for any energy de- 
pendence, and the results are identical with those shown in figure 2. Consid- 
ering the strong similarity between our uncorrected results and those of Brink- 
man et <sZ., the fact that the data were taken with the same instrument, and 
the fewer number of observations, which would lead to greater uncertainty in 
the correction factors, we suspect that the latter experiment is consistent 
with no energy dependence. 


III. CROSS-CORRELATION ANALYSIS 


The cross correlation between the counting rate S in the low-energy 

J_i 

band (2. 1—5. 2 keV) and the counting rate S in the high-energy band (5.2—16.4 

K 

keV) was also examined. For each observation of Cyg X-l the cross correlation 
function (CCF) 


CC(u) = 


N - u 

l 

i = 1 


< S L,i - V < S 


H , i + u 


-vi 


q L,i °H,i + u 


[ r L (°) r H (°)] 


1/2 


( 8 ) 


. *• . a T . a„ . 

1 = 1 L,1 H, l 


was calculated; this function is normalized to fall between -1 and +1. As be- 
fore, the CCF will only represent the behavior of the source in the limit where 
the fluctuations due to counting statistics are negligible. By an argument 


analogous to that leading to equation (6) we find 


CC(u) = 


V V 
L,x H,x 


V V 
L H 


1/2 


CC (u) 
x 


(9) 
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where CC is the true CCF of the source with variances V and V in the 
* L,x H,x 

low- and high-energy bands, and V L and V are the detected total variances in 
these bands. In practice, we calculated an average CC(u) for the 71 observa- 
tions and then corrected it by the square root of the ratio of the average 
vaiiances to obtain CC^(u) . This corrected CCF is shown in figure 4 . 

ine results of the cross-correlation analysis shown in figure 4 indicate 
effectively a 100 percent correlation in time between the high- and low-energy 
data (at least on time scales ~ 0.2 s) with no indication of any asymmetry. 
That is, the high-energy events neither lead nor lag the low-energy events. 

The most natural interpretation of this result is that the underlying shots, 
hot spots, or flares that give rise to the high- and low-energy photons are 
identical and have essentially the same spectrum from flare to flare. (If the 
photons In these two bands arose from different regions of the source, or the 
spectrum of a hot spot was quite variable, then any significant cross correla- 
tion would be unlikely.) A further argument that the sources of the low- and 

high-energy photons are indeed identical is that the ACF and CCF are identical 
in shape. 

IV. CONCLUSIONS 

We have shown that for data obtained in the time interval 1972 January 
9-27 the short-term time variability of the source Cyg X-l can be represented 
by a shot-noise model with characteristic time constant of 0.45 s, that the 
time variability is independent of energy, and that the cross correlation be- 
tween events in different energy bandwidths is as large as 100 percent for 
time scales ~0.2 s. These results are consistent with the production of all 
of tne 2 16— keV X— rays in random flares (of lifetime —0.45 s) in the accretion 
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disk around the compact object. Finally, we have shown that data obtained in 
a previous observation of the source (Brinkman et at. 1974) in 1971 June are 
entirely consistent with these results, indicating possibly that the short- 
term time variability seems to be constant on time scales on the order of at 
least several months. This conclusion is further substantiated by a recent 
rocket observation of Cyg X-l in 1973 October (Boldt et cd . 1974). 

We note that all of these observations were performed after the spectral 
transition of 1971 March. A comparison of the variability before and after the 
spectral transition may provide further insight into the emission mechanism 
of this object. 
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tory. This work was supported by NASA under grant NGR 33-008-194 . This 
Letter is Columbia Astrophysics Laboratory Contribution No. 107. 



10 


REFERENCES 

Boldt, E. , Holt, S., Rothschild, R. , and Serlemitsos, P. 1974, "What Is 

Special about Cygnus X-l," NASA Goddard Space Flight Center, Green- 
belt, Md., preprint X-661-74-323 . 

Brinkman, A. C. , Parsignault, D. R., Schreier, E., Gursky, H. , Kellogg, E. M. , 
Tananbaum, H., and Giacconi, R. 1974, Ap. J. , 188 , 603. 

Giacconi, R. , Kellogg, E. M., Gorenstein, P., Gursky, H. , and Tananbaum, H. 
1971, Ap. J. (Letters) , 165, L27. 

Jenkins, G. M. , and Watts, D. G. 1968, Spectral Analysts and Its Applications 
(San Francisco: Holden-Day, Inc.), Ch. 5. 

Terrell, N. J., Jr. 1972, Ap. J. (Letters) , 174, L35. 

Thorne, K. S . , and Price, R. H. 1975, Ap. J . (Letters) , 195 , L101. 



IX 


FIGURE CAPTIONS 

FIG. 1 . — a) Data obtained by the l/kuipu satellite in a single obser- 
vation of Cyg X-l on 1972 January 9. b) The same data with the triangular 
response of the collimator removed. The error bars shown are based on count- 
ing statistics. The portion of the data indicated by A — A was selected for 
further analysis. 

FIG. 2. The average of 71 autocorrelation functions and power density 
spectra. Only the first 36 values of the autocorrelation function are shown. 
The error bars are based on the scatter of each data point about the mean for 
the results uncorrected for the bias produced by counting statistics. The 
error bars shown for the corrected functions also include the uncertainty in 
the correction factor. 

FIG. 3. -The corrected average autocorrelation functions for data in 
two energy bandwidths. Only the first 36 values of the autocorrelation func- 
tions are shown. The error bars shown are those for the high-energy data only. 
The error bars for the low-energy data are only slightly smaller. The error 
bars are based on the scatter of the individual data points about the mean 
and the uncertainty in the correction factor. 

FIG. 4. — The average cross-correlation function corrected for bias 
produced by counting statistics. The error bars are based on the scatter of 
the individual data points about the mean and the uncertainty in the correction 
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